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We provide an exact finite temperature extension to the recently developed Riemann-Hilbert 
approach for the calculation of response functions in nonadiabatically perturbed (multi-channel) 
Fermi gases. We give a precise definition of the finite temperature Riemann-Hilbert problem and 
show that it is equivalent to a zero temperature problem. Using this equivalence, we discuss the 
solution of the nonequilibrium Fermi-edge singularity problem at finite temperatures. 



I. INTRODUCTION 

The response of a Fermi gas to an external, time- 
dependent perturbation is a well-known hard problem. It 
has repeatedly been studied in various contexts, mostly 
on systems exhibiting the Orthogonality Catastrophe, 
Fermi-edge singularities, or the Kondo effecti The 
technical difficulties are connected to the necessity of 
summing an infinite number of logarithmically diver- 
gent diagrams, and has led to the development of a 
multitude of methods to tackle this task, including dia- 
gram summations, singular integral equation approaches, 
bosonization, and summation schemes based on Slater 
determinants^. Many of these approaches are technically 
involved or require assumptions such as, for instance, 
the separability of the external potential. In addition, 
the recent controversies concerning nonequilibrium Fermi 
systems^ have shown that extensions of these methods 
are problematic and hard to achieve. A unified approach, 
depending on a minimum of assumptions is, therefore, 
desirable. 

Very recently, a major step in this direction has been 
made in the form of the Riemann-Hilbert approach£&, 
and applied to the calculation of Fermi-edge singular- 
ities at T = in Fermi gases under nonequilibrium 
constraints^. In this approach, the calculation of gen- 
eral response functions in the perturbed Fermi gas is 
transformed into the problem of solving a (generally non- 
abelian) Riemann-Hilbert (RH) boundary value problem: 
The response functions are rewritten in terms of complex 
matrix functions Y(t) that satisfy boundary conditions, 
for Imt — * ±0, expressed by a matrix function R(t). In 
most cases of interest, R(t) is directly related to the scat- 
tering matrix of the time-dependent potential. 

In this paper, we extend the RH technique to finite 
temperatures. We first give a precise definition of the fi- 
nite temperature RH problem. We then show that the fi- 
nite temperature problem is equivalent to a different zero 
temperature RH problem through the bijective mapping 
of the time variable t i— » tanh(5t), with (3 = l/fc^T. In 

simple (i.e. abelian) cases, this allows us to construct the 
explicit solution, but generally the resulting RH prob- 
lem remains highly nontrivial. We show, however, on 
the co ncret e example of the nonequilibrium problem of 
Refs. I3 EI5I how the equivalence between finite and zero 



temperature solutions can be explored to construct ap- 
proximate solutions to the corresponding RH problem in 
the high and low temperature limits. 

The paper is organized as follows: In Sec. ^ we sum- 
marize the essential steps in the zero temperature RH 
approach. In Sec. lIIII the finite temperature RH problem 
is defined, and the equivalence to the zero temperature 
case is shown. Finally, in Sec. IIVI we discuss the finite 
temperature extensions for some simpler examples and 
for the nonequilibrium Fermi-edge singularity problem 
discussed in Refs. l3l lllril 

II. ZERO TEMPERATURE CASE 

We start with the introduction of the required concepts 
and notations for the zero temperature RH approach, 
following Ref. 0. 

The approach has been developed for the calculation of 
response functions in a nonadiabatically perturbed sys- 
tem of fermions. We consider a gas of noninteracting 
fcrmions that is exposed to a time-dependent potential 
M. The Hamiltonian is assumed to take the form 

H(t) =H +J24M(t,e,e')a e ,, 

(1) 

H = ^ea^a £ , 

e 

where e labels the single-particle energies of the electrons, 
and a| is a vector of the electron operators a\ ■ , where 
j = l,...,N is a channel index labeling any further 
classification of the single particle states. The potential 
M(t, e, e) is a N x N matrix in channel space, which we 
assume to be nonzero only for some finite time interval 
< t < tf, during which the external perturbation is 
switched on. 

The quantities of interest are response functions of the 
type 

XR={R). (2) 

At zero temperature, (...) denotes the average over the 
ground state |0); at finite temperatures, as considered 
in Sec. IIIII it represents the thermal average over the 
noninteracting system states at t = 0. The operator 
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R is related to the time-evolution operators U(tf) and 

Uo(tf) of Hit) and H , respectively, or equivalently to 
the scattering matrix S for the potential M. For the 
present discussion, the precise form of R is unimportant. 
A summary of necessary conditions on R is given in Ap- 
pendix |A] while a detailed discussion can be found in 
Refs. |!|0and|(| some examples are given at the end of 
this paper. 

Since the Hamiltonian is quadratic in the electron op- 
erators, the action of the many-body operator R is fully 
specified by its action on single-particle states (see the 
appendices). Let R = {Rj,j>(e, e')} De the matrix formed 
of the matrix elements of R between the single-particle 
states (j, e), Rj.j>{e,e') — (\a e jRa e >j>\), with |) the true 
vacuum. For the following treatment, it is essential to 
assume that R is diagonal in the time representation, 
R = R(t) (with t the variable conjugated to the energy 
representation e) . For the response functions relevant for 
Fermi-edge singularity problems, it can be shown&iS (see 
Appendix that this assumption is valid as long as the 
variation of the scattering matrix is slow compared to the 
delay time of the scattering process itself. 

For the noninteracting system, |0) is a Slater deter- 
minant in the single particle states. As shown in Ap- 
pendix [51 we can write 



X R = det((l-f) + fR), 



(3) 



where R represents the matrix elements of R taken be- 
tween the single-particle states, and / is the Fermi func- 
tion, which has the matrix elements 



(4) 



In this expression, we have set the chemical potential 
/i to zero. This is possible without loss of generality 
even in the nonequilibrium case of several Fermi seas 
j with different chemical potentials fij, such as for bi- 
ased tunneling barriers. For the unperturbed system 
with uncoupled Fermi seas, we can shift the single- 
particle energies by e -* e — [ij and set all chemical 
potentials to zero. In other words, the fermion opera- 



tors are replaced by a e 



e w *a e 



The Hamiltonian 



keeps the same form except that the potential M(t), and 
thus R(t), acquires an additional dependence on time 3 , 

The logarithm of the determinant in Eq. © can be 
written in the form 



lnxfl=Tr[/ln#] 



In 



-Tr[ln(l -/ + /#)-/ In fl] 

(2) (5) 
Xr, 



with Tr being the trace over channel indices j and en- 
ergy. ln%^ = Tr[/lni?] consists of the diagonal terms 
of lni? summed over the occupied states of the unper- 
turbed system. Depending on the system under investi- 
gation, it expresses, for instance, a shift of the ground 
state or threshold energy (Fumi's theorem^ if R is the 



scattering matrix) together with an imaginary correction 
for nonequilibrium conditions^, or the average transfer 
of charge across a barrier 3 . The second term, ln%^ , 
contains the nontrivial effects close to the Fermi surface 
which arc, for instance, associated with Fermi-edge sin- 

(2) 

gularities. In the following, we will focus on Iv-Xr only. 

It is convenient to switch from the energy represen- 
tation to the time representation of the trace, in which 
R is diagonal. The Fermi function, however, becomes 
nondiagonal in time and reads 



(6) 



(Here and henceforth we set h = 1.) The logarithm of 
an infinite matrix is represented by introducing a depen- 
dence on a parameter A as 



R — > R(X) = exp(Alni?) 



(7) 



and by writing 
-l 



1 (2) 



dA / ditr 



[(l-f + fRy'I-fR- 



dR 
dA 

(8) 

where tr denotes the trace over the channel index j only. 

The difficult task is to find the inverse of (1 — / + fR). 
For this, a complex N x N matrix function Y(t) is in- 
troduced, which solves an auxiliary RH problem, defined 
by: 

(a) Y(t) is a nonsingular N x N matrix function, which 
is analytic and nonzero in the complex plane, except 
on the interval (0,tf), where it satisfies the boundary 
condition 



Y_(t)Yz 1 (t) = R(t), 



(9) 



for Y±(t) — Y(t ± iO), and except at the end points 
t e = 0,tf, where Y(t) may vanish or be weakly sin- 
gular as \Y(t)\ ~ \t-t e \ £ for -1 < e < l,e ^ 0. 



(b) Y(t) has the asymptotic condition 
Y(t) -> 1 when |i| -> 



(10) 



From these two properties, two important identities fol- 
low: 



fY+f = Y+f, 
/y_/ = fY_, 



(11) 
(12) 



which equally hold if Y is replaced by Y 1 . In the time 
representation, Eq. (jl 1|) reads 



i/2n , s i/2n , , i/2t: 



h-t + iO w t-t 2 + «0 



h -t 2 + i0' 

which is proved by closing the contour in the upper half- 
plane. A similar contour in the lower half-plane proves 
Eq. O- 
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In terms of the functions Y±, the inverse of (1— f+fR) 
is expressed by 

(i - / + fny 1 = y+ [(i - 1 + fYi 1 ] , (w) 

which can be verified by multiplication by (1 — / + fR) 
from the right or left, and by making use of the identities 
(fTTfl and (15)1. 
Finally, the following result 



dttr{[A, f]B} 



dt lim tr — 

t'-kt 27T 



A(t)-A(0 



t - f + iO 



B{t') 
(15) 



for any differentiable matrix functions A and B, allows 
to write Eq. JHJ) in the compact form 



1 (2) 



, dA / dt tr . — — , 
2ir J J \ dt + dA 



dYl 



-YZ L R 



- r dR 



(16) 



The response function Xr i therefore, is fully expressed 
in terms of the solution Y(t) of the RH problem. 



III. FINITE TEMPERATURE CASE 

We show in this section, that we can define a finite 
temperature version of the RH problem such that the 
important formula i|14|) (and, therefore, Eq. (jlGf! ) remains 
valid. Subsequently, we show that this new RH problem 
is equivalent to a zero temperature RH problem. 



A. Finite temperature RH problem 

Note that Eq. (|14fl remains valid if we generalize the 
identities (|TT]) and (T5)l to 



fY + f = Y + f + F 1 

fY.f = /y_ + f, 



(17) 
(18) 



with the same F — F(t,t'), a (finite) matrix function, 
appearing in both identities. 

At finite temperatures T, the Fermi function becomes 



The Fermi function is antiperiodic in the imaginary 
time direction, f(t + i(3,t') = —f(t,t'). If we impose 
that Eq. and JTBJ are invariant under such trans- 
formations, Y(t) and F(t,t') must be (anti)periodic in 
the imaginary time direction as Y(t + i(3) = Y(t) and 
F(t + 0,t') = F(t,t' + 0) = -F(t,t'). This allows us to 
restrict the analysis to the strip —(3/2 < Imt < (3/2. 

The finite temperature RH problem can then be formu- 
lated as follows: 

(i) Y(t) is a nonsingular N x N matrix function, which 
is analytic in the strip Sp = {t € C| — (3/2 < 
Imt < (3/2}, except on the interval (0, tf), where 
it satisfies Eq. @, and except at the end points 
t e = 0,tf, where it may vanish or be weakly singu- 
lar as \Y(t)\ ~ \t - t e \ £ with -1 < e < 1, e ^ 0. 

(ii) Y(t) tends to definite finite values as \t\ — > oo in Sp 
(the values on the far right or far left, however, are 
generally different). 

(iii) Y(t) can be analytically continued through t±if3/ir, 
for real t, such that Y(t + i(3/2) = Y(t - 0/2). 
For convenience, we normalize Y(t) (by multiplica- 
tion by a constant matrix from the right) such that 
Y(±0/2) = 1. 

Condition l)ui)l is absent in the zero temperature case, but 
is essential for the validity of Eqs. (|17fl and l|18|) . 

Indeed, let us first focus on Eq. I|17|) . We integrate 
fY+f over the contour C = C\ + ■ ■ - + Cs shown in Fig.^ 
In the figure, r is a large number, which we eventually 
send to infinity. Since Y(t) remains bounded as \t\ — > oo, 
the integrals over the vertical lines C2 and C5 vanish in 
this limit, and we obtain Eq. (|17|) with F given by the 
integral of fY+f over the line through i(3/2, parallel to 
the real axis. 

A similar contour in the lower part of Sp leads to 
Eq. (18)1 . with F given by the integral over the line 
through —i/3/2, parallel to the real axis. Both integrals 
for F are identical due to the periodicity and analyticity 
of Y(t), and Eqs. (33 and (HHJ) are fulfilled. 

The solution of the finite temperature RH problem, 
together with Eqs. i|17|l and l|18|) . immediately leads to 
the solution for the response function \R by means of 
Eq. i|ltj|) . (Note that Eq. I|15|l remains true when replac- 
ing (t - t')- 1 by [(jS/tt) sinh(^ - t')/P)]~K) 



f{t,t') 



i/2n 



f si nh(|(t-t')+i0) 



(19) 



with (3 = l/k B T. With this function, Eq. JUJ reads 
explicitly 



dt 



i/2/3 



-Y+(t) t- 



i/2/3 



sinh(|(ti -t + iO)) +v sinh(2(t - t 2 + iO)) 
i/2/3 



sinh(5(ii -t 2 +*0)) 



+ F(t u t 2 ). (20) 



B. Equivalence to zero temperature case 

The finite temperature RH problem, defined by con- 
ditions (fif) — (fTTT|> . is equivalent to the RH problem at zero 
temperature, defined by conditions @ and (0, upon a 
redefinition of the boundary condition © . 

The equivalence is shown by means of the mapping 



1 1 — ► vd — tanh(-|i), 



(21) 



which is often used to switch between zero and finite 
temperatures, and which describes a bijection between 
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the strip Sp and the slit plane S = {w S C\w $ 
(— oo, — l],w & [l,oo)}. Let us set Y(w) = Y(t) for 
w = tanh(-|t). 

• The boundary condition © for Y(t) and t <E (0,i/) 
becomes a boundary condition of the same type 
for Y(w) and w € (0,w/), < w/ = tanh(5t/) < 
1, and the matrix R(w) = R(t). 

• The analyticity of Y(t) in Sp\[0,tf] ensures the 
analyticity of Y(w) in 5\[0,io/]. The periodicity 
Y(t + 0/2) = Y(t - 0/2) for real i is nothing but 
the condition of analyticity of Y(w) through the 
real axis at \w\ > 1. 

• t = ±0/2 maps onto infinity in the w— plane. With 
the normalization Y(±if3/2) = 1, we have Y(w) — » 
1 as \w\ — » 00. On the other hand, the boundedness 
of Y(t) as \t\ — ► 00 in Sp implies that Y(w) is finite 
at w — ±1. 

This proves that the zero and finite temperature versions 
of the RH problem are equivalent, upon the change of 
boundary condition R(t) — > = 

With the mapping, the identities (|17|l and (|18|l are 
mapped onto Ijlll) and 112fl . too. The infinite ^-integrals 
over the real axis, map onto (—1,1) in the w-plane (in- 
cluding the cut [0, Wf]), while the integrals over the con- 
tours C3 and C4 provide the remaining parts of the real 
axis in the w-plane, (1, 00) and (— 00,— 1), respectively, 
at which Y(w) is analytic. The integrands can be iden- 
tified through the formula 

sinh(|(* - i')) = cosh(|t) cosh(|i')(w - w'). (22) 

The dependence on cosh(Si) in the integrand cancels 
with ^ = (7r//3)/cosh 2 (|t). After factorization of 
(7t//3) cosh(-gti) cosh(gt 2 ), we see that Eqs. (JT71 and JTSJl 
are identical to (fTTft and ( fT5f> . 



IV. EXAMPLES 

A solution to the finite temperature RH problem allows 
the calculation of the response function (|16|l . We illus- 
trate this with a few simpler examples of how the equiva- 
lence to the zero temperature RH problem can be used to 
construct explicit finite temperature solutions. We then 
discuss approximate finite temperature solutions for the 
nonequilibrium RH problem of Ref. 0. 



A. Some simpler examples 

The RH problem appeared early in the study of Fermi- 
edge singularities. It is part of the solution of the Dyson 
equation for the Green's function, formulated by Nozieres 
and De Dominicis as a singular integral equation^, for 



which the standard method of solution consists in the 
transformation into a RH problem 9 . In Nozieres and 
De Dominicis' case, the boundary condition R is a con- 
stant, for t E (0,tf), given by the single-particle single- 
channel scattering matrix R — S — e~ 2lS , where S is the 
scattering phase shift (\S\ < it/2). The solution is given 
by 



Y(t) = ((t-t f )/t) 



S/-K 



(23) 



The finite temperature extension to this expression was 
achieved by Yuval and Anderson^ in their application 
of the Nozieres-De Dominicis approach to the Kondo 
problem. It resulted from writing l/sinh(7r£//3), the free 
Green's function in the kernel of the integral equation, as 
the sum of (— l) n /(f + in(3) over integer n in the kernel 
of the integral equation. The corresponding RH problem 
becomes the product of zero-tempcrature RH problems 
with branch cuts shifted by in(3. Hence, the finite tem- 
perature solution found by Yuval and Anderson reads 



nt) - n 



t + in/3 



8/7, 



sinh(ft) 



(24) 

With the replacement t 1— > tanh(^i) in Eq. (|23|l and the 
use of Eq. I|22|l . this solution can immediately be veri- 
fied, as constant matrices R and R coincide. Note that 
with the mapping we obtain an additional constant fac- 
tor cosh(^i/)' 5 / 7r , reflecting the fact that a solution of the 
RH problem remains a solution upon multiplication by 
a constant matrix from the right. With our normaliza- 
tion ijmjl this factor is suppressed. Such constants are 
unimportant since they drop out in the physical solution 

Eqs. Il2.'j[| and (|24() remain valid if R is a time- 
independent diagonalizable matrix. In this case, 5 be- 
comes a matrix with eigenvalues representing the phase 
shifts of the multi-channel scattering matrix R. 

For scalar, but time-dependent R(t), the solution of 
the RH problem can be constructed in a standard way 9 , 
and is given by 



n*)=ex P( - 



t'-t 



(25) 



At finite temperatures, the mapping allows here the ex- 
plicit solution, too, which has the same form as the pre- 
vious equation, with R replaced by R, 

( 1 /• tanh (? 4 /) ]nR(w') \ 

Y(t) = exp / dw'— \ ' r 

w y \2TTiJ Q w' - taah(%t)J 

1 f*f dt , cos Hp) hiR(t') \ 



exp 



2f3iJ cosh(|t')sinh(|(i'-<)) 



(26) 

Notice that all these expressions become singular at t — > 
or t — > tf. These divergences are spurious and due to 
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the approximation of an infinite bandwidth (given by the 
form of the Fermi function) , and have a natural cutoff in 
the true inverse bandwidth 8 , l/£. 

In these simple examples, the response function <|16|) 
can explicitly be calculated. Let us, for instance, consider 
Eq. (22J. Since S _1 4§ = -2i6, we obtain 



^X<k\tf) = 



tr{5 2 } 
2tt 2 



f , d , 
dt— In 

di 



sinh(S(t-i/ + iO)) 



tr{(5 2 } 



In 



^ 

sinh(|t + iO) 

sinh(7rij//3) 



sinh(7r//30 



(27) 



again with the cutoff at the inverse bandwidth l/£. For 
short times tf <C (3, we obtain a power-law similar to the 



well known zero temperature result Xr (tf) 
which, however, crosses over into an exponential decay 
as Xn(t f )^ e -( s2 ^m. 



1. Zero temperature case 

In the nonequilibrium Fermi-edge singularity 
problem 4,5 , R(t) is a unitary time-dependent scat- 
tering matrix coupling two Fermi liquids, j = 1 and 
j = 2, at the left and at the right of a tunneling 
barrier. The nonequilibrium constraint is introduced 
by a finite voltage bias V, keeping the left and right 
chemical potentials at the values [i\ = V + [i and 
fj,2 = \i, respectively (the electronic charge e is set to 
unity). As noted in Sec. [H] we can shift /j,± t 2 to in 
the unperturbed system by passing to the interaction 
representation Hq — ^ Hq = Hq — fiiNi — where 
are the number operators for the electrons in 
the electrodes 1 and 2. This is equivalent to a gauge 
transformation, transforming the electron operators as 



'< ,1 



= e lfij a. 



The voltage difference now ap- 
pears only as an additional time dependent phase factor 
in the matrix elements of R(t) that couple electrodes 1 
and 2: 



B. Biased tunneling barriers 

As a more involved example, let us consider the non- 
abelian RH problem, in which R(t) is a time-dependent 
matrix that does not allow a time-independent diagonal- 
ization. Such situations naturally arise in nonequilib- 
rium systems, and have recently been studied with the 
RH approach 4 * 5 -: Consider, for instance, the coupling be- 
tween a defect state in a tunneling barrier and the con- 
necting electrodes that are kept at a voltage bias V. The 
defect is assumed to be either in its ground state g or 
in in an excited state e. We further assume that in the 
ground state the defect is fully screened and does not 
interact with the conduction electrons, whereas in the 
excited state it acts on the conduction electrons through 
the potential M. The diagonal matrix elements of this 
potential describe the scattering of the electrons in the 
electrodes, while the nondiagonal matrix elements rep- 
resent a modification of the tunneling rate. The matrix 
R = S is the scattering matrix for M. The response 
function Xn(t — tf ~~ 0) 1S related to the transition rate 
of the defect into the excited state, p g ^ e (AE), by 



dt xn{t) e 



iAEt 



(28) 



which is the Golden Rule expression rewritten as a time 
integral. AE is the bare energy difference between 
ground state and excited state of the defect, and u is 
the bare amplitude in the Hamiltonian for the transition 
g -> e. 

A closed solution to the resulting time-dependent RH 
problem is not known. In the short time limit t -C 1/V 
and in the long time limit t 1/V, however, it is possible 
to construct approximate solutions^* 5 ". We first give an 
example for such an approximation at zero temperature, 
and then discuss how it extends to finite temperatures. 



R(t) 



R11 e lVt R 12 
~ lVt R2\ R22 



(29) 



where i?n and R22 are the backscattering matrix ele- 
ments, and R12 and R21 describe transmission through 
the barrier. For V = 0, the scattering matrix R is time- 
independent during (0, tf), and the explicit solution of 
the RH problem is given by Eq. I|23fl . with S the matrix 
of phase shifts, R — exp(—2iS). 

If V 7^ 0, an explicit solution to the RH problem is 
not known. The naive solution, obtained by multiplica- 
tion of Eq. (|23|1 by ^-dependent phase factors exp(±iVt) 
would solve Eq. 10 but violate condition 0, expressed 
by Eq. I|1(J|I . because of the divergence of exp(±iVt) for 
large complex-valued arguments t. Such a divergence 
cannot be compensated by multiplication by t-dependent 
analytic matrix functions from the right. The solution 
would become invalid. An approximate solution, 
however, can be constructed as follows^ 4 * 5 -: For V > 0, 
we decompose R(t) into the product 



R(t) 



1 

a/z 1 



a 
b 



1 a z 
1 



(30) 



with z = exp(iVt), a — i?n, b — det R/Rn = 1/R 22 , 
a = R21/R11, and a' = -R12/-R11. Now define the func- 
tion 



,'lnD, 
ip(t) = exp — — In 
1 2m 



t 



t - t 



(31) 



where D = diag(a,6). This function satisfies the bound- 
ary condition 

^.(t)i;- 1 (t)=D (32) 
at t € (0, tf). To approximate the unknown solution 
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Y(t), we consider the following function W(t): 



2. Low temperatures: ksT -C V 



w 1 (t)=m, 



W 2 (t) = 



1 — a'z 
1 



W 3 (t) 



1 

a/z 1 



tf>{t), W*(t)=1>(t), 



where the functions W s are defined in the regions s — 
1, . . . , 4 indicated in Fig. |U (a). W has branch cuts be- 
tween the boundaries of all 4 regions s. Most importantly, 
between regions s = 2 and s = 3, it satisfies 



w 3 (t)w 2 -\t) 



W-(t)W7 1 (t)=R{t). (34) 



Furthermore, W(t) — > 1 as \t\ — > oo. Hence W(i) satisfies 
the conditions © and (jbj) of the RH problem, at the price 
of the additional vertical branch cuts through and t / . 
These branch cuts, however, become exponentially small, 
oc e~ v ^\, with increasing distance to the real time axis. 
For instance, consider the vertical cut at Ret — in the 
upper half-plane. We have 



W 2 (t)W^(t) 



1 



1 -a'e-^l 
f 

+ 0(e- v '*l). 



(35) 



Similar exponential suppressions hold for the remaining 
vertical cuts. Therefore, Y(t) = Wif) is a good approxi- 
mation to the solution of the RH problem if we consider 
times t such that \t\, \t — tt\ ^> 1/V. Notice that in the 
boundary condition i|32|) only a — R%i and b = 1/RW 
are involved, the matrix elements of backscattering into 
the left electrode at fix = V + fi and into the right elec- 
trode at /i2 = /i. This reduction of the full scattering 
matrix can be explained in a simple picture based on the 
uncertainty relation: For long times t 1/V, the en- 
ergy uncertainty Ae of the scattered particles is smaller 
than the voltage, Ae <C V. The Fermi surface of the 
other electrode, therefore, is not seen by the scattered 
particle, and only the backscattering part of the scatter- 
ing matrix is of importance. For the left electrode, this 
is Rn, while for the right electrode, the Pauli principle 
forbids virtual trajectories into the left electrode. Af- 
ter removal of these, the scattering amplitude becomes 

R22 — R21R11R12 = 

Since Rxi and R22 alone are not unitary the real phase 
shifts S are replaced by complex numbers. The implica- 
tions, mainly, a broadening of the Fermi-edge singularity, 
are s hort ly discussed below. For more details we refer to 
Refs. 

In the limit \t\ or \t — tf \ <C 1/V, on the other hand, 
nonequilibrium constraints play no major role, and the 
equilibrium expression (|23|l is expected to become appro- 
priate. This indeed is confirmed by the same argument of 
the uncertainty relation: For small times, t <C 1/V, the 
energy uncertainty Ae of the scattered particles becomes 
large, Ae 3> V, and both Fermi surfaces at fi\ and fi2 are 
coupled. The full scattering matrix is involved. 



Let us focus now on low temperatures, fc^T <C V. 
An approximate solution to Y(t) to the RH problem can 
be constructed in the same way as for the zero tempera- 
ture case: We decompose R(t) according to Eq. PU|) . and 
consider the RH problem for ip(t), described by Eq. 132fl . 
Since this problem is diagonal with a time-independent 
boundary condition, it can be exactly solved by the map- 
ping on the zero-temperature problem. Its solution is 
given by Eq. I|24|) with — 5/ir replaced by ln_D/2i7r, 



W) = CX P I ^— In 



sinh(§£) 



(36) 



The function W(t) is given by the same conditions l|33(l as 
before, but where the regions s = 1, . . . , 4 now are defined 
in the strip Sp, with additional branch cuts at Imf = 
±/3/2 (see Fig. 121(b)) The validity of the approximation 
on the vertical branch cuts remains unchanged from the 
zero temperature case, and corresponds to \t\, \t — tf \ ^> 
1/V. On the additional horizontal cut, W2(t + i(3/2) and 
W 3 (t — i(3/2) differ in their nondiagonal parts, which are 
of the order of e~@ v / 2 . This leads to corrections for the 
identities (|II|I and (|12|) of the form 



fW + f = W+f + F + c-^ 2 F + , 
fW-f = fW_ +F Q + c-W 2 F_, 



(37) 



where Fo is the integral over all coinciding parts of fWf 
along Imt = ±/3/2. F + 7^ F_ are the contributions due 
to the mismatch of W% and W3. Both, F+ and F-, are 
finite quantities of the order of unity. Since (3 ^> 1/V, 
their contribution is negligible compared with the error 
arising from the vertical cuts. 

On the other hand, for |i| < 1/V < /3, the volt- 

age (i.e. the nonequilibrium situation) as well as the tem- 
perature become unimportant. The solution to the RH 
problem, therefore, is that of a zero temperature equilib- 
rium system, expressed by Eq. (|23|l . 



3. High temperatures: ksT ~S> V 

At high temperatures, fc^T ^> V, we expect that V 
plays no longer an important role. Indeed, since (3V <C 1, 
the naive solution of the RH problem, 



Y(t) = 



sinh(f^) 



smh(%(t-t f )) 



(38) 



now becomes accurate. Because of the restriction of t to 
the strip Sp, Y(t) cannot diverge exponentially, so that 
\Y(t)\ remains bounded as \t\ — > 00 in Sp. Yet, Eq. 
violates the condition of periodicity, Y(t + i(3/2) ^ Y(t — 
i(3/2) due to the exponentials e ±lVt . Again, this leads to 
different integrals F from the integration along t + i(3/2 
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and t — respectively. Since (3V <C 1, however, we 

can expand the exponential functions, and we are led to 
similar identities as before 



fY + f = Y + f + F +pVF+, 
fY_f = fY_+F +pVF_, 



(39) 



with F± of the order of unity. The solution Y, and hence 
the inverse of the matrix, which is solved by the RH 
problem, as well as the logarithm of the response function 
are accurate within corrections of the order of (3V . 



4- Nonequilibrium Fermi-edge singularities 

The results above allow us to extend the nonequilib- 
rium Fermi-edge singularity problem to finite tempera- 
tures. We consider the case where R(t) has the form 
given in Eq. flU with R = S = e~ 2lS . 

The high temperature limit leads to an expression, 
which is unchanged from the equilibrium result l|27|l. 
At lower temperatures fc^T <C V, on the other hand, 
we can use the approximation Y(t) w W(t) for for 
\t\, \t — tf\ ^S> 1/V. With these functions and the decom- 
position of R in Eq. H3()|) . the matrix products in Eq. I |16(l 
can be evaluated to 



tr 



dt + 



_id_R 
dA 



= 9(t) 



d tr{(lnL>) 2 } 



dA 



2m 



where 



9(t) 



dt 



hi 



sinh(|(t + iO)) 
nhm-t f +iO)) 



-iVC(X), 
(40) 

(41) 



and C(A) is a time- independent function of A. The term 

( 2) 

proportional to C(A) leads to a contribution to Inxk 
that is linear in Vtf, and hence consists in a voltage 
dependent (generally complex) shift of the ground state 
energy, similar to that of ln%^ . The time integration 
over g(t) has to be cut off at t — 1/V and tf — t = 1/V, 
where the solution Y(t) crosses over into the equilibrium 
solution. The remaining integral can be estimated by 
integrating the equilibrium solution from t = l/£ to 1/V. 
This leads to 



(2») 2 



3inh(7rf / / 0) 



sinh(7r//3V) 



tr{<5 2 } 



In 



sinh(7r//3V) 
sinh(7r/^) 



(42) 



where we recall that tf > 1/V. Since 1/f < 1/V -C /?, 
this can be rewritten as 



ln X g)(t / )=- % V C 7-M^ln 



5mh(irtf //3) 



(2tt) 2 [ 7T//3V 

-^>ln«/V). (43) 



For l/£ <C tf <C 1/V, this function crosses over into the 
equilibrium solution lt77|) . 

Since both S and hiD are time-independent (for t S 
(0,tf)) we see that, in this limit /cgT <C V, the exten- 
sion from zero to finite temperatures, originally achieved 
by the mapping t i— > tanh(7rf//3), can here be obtained 
by a rule of thumb of replacing tf by (/3/n) sinh(7rf///3) 
in the logarithms. For XR.{t) this means that the zero 
temperature power laws of the form tj are replaced by 
[smh(7rf///?)] Q . 

Therefore, finite temperatures replace the power-law 
tails by an exponential decay with a decay time propor- 
tional to the inverse temperature (3. In other words, while 
the system at zero temperature keeps a long-term mem- 
ory to the initial perturbation, finite temperature fluc- 
tuations suppress it within a characteristic time set by 
0/a. 

The transition rate p g ^ e (AE) is determined by the 
Fourier integral (|28p. An explicit calculation of the in- 
tegral is possible and leads to a Beta function, yet the 
following estimates are physically more instructive: 

At zero temperature both, the singularity at t — > as 
well as the power-law tails for t — > oo contribute equally 
to the integral. In the equilibrium situation, we obtain 
the well-known Fermi-edge singularity, a power-law de- 
pendence on energy, p g ^ e ~ AE 4 '" "* _1 . For V ^ 0, 
AE and the exponent, acquire an imaginary correction. 
The nonequilibrium effects dominate for small energies 
AE < V, and the imaginary corrections to AE lead to 
a broadening of the Fermi-edge singularity. The high 
energy tails, AE ^> V, coincide with the equilibrium so- 
lution. 

At finite but small temperatures, fc^T <C V, the ex- 
ponential decay of the response functions for t ^> l/fc^T 
causes an additional, Lorentzian, broadening of the sin- 
gularity for AE < ksT. At high temperatures, ksT ^> 
V, the nonequilibrium effects are completely washed out. 
For small energies, AE <C fc^T, the exponential decay 
then dominates the integral, and the singularity is broad- 
ened to a Lorentzian with a width proportional to fc^T. 
On the other hand, for high energies AE ^S> ksT, the 
exponential decay is slow compared to the quickly oscil- 
lating phase factor. The behavior of the integral can be 
estimated through a saddle point approximation, leading 
to p g ^ e (AE) ~ exp(— AE/ksT), corresponding to an 
activated transition such as for a noninteracting system. 

In summary, nonequilibrium constraints and tempera- 
ture have the similar effect of a broadening of the Fermi- 
edge singularity with two different shapes: A Lorentzian 
broadening from the temperature with a width given by 
ksT, and a Lorentzian to the power of an expression 
depending on In D from the nonequilibrium constraints 
with a width determined, to first order, by V. Away from 
the resonance, but still for AE <C fc^T, the transition 
rate has a power-law dependence on energy, which even- 
tually turns into an activated behavior for AE > ksT. 
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V. CONCLUSIONS 

We have shown that the RH approach developed in 
Refs.iandi has an exact finite temperature extension. 
This finite temperature case can be mapped onto a differ- 
ent zero temperature problem, for which much is math- 
ematically knownii, and which allows, in some simpler 
cases, an exact solution. Based on this equivalence we 
have constructed approximate solutions for the nonequi- 
librium RH problem at finite temperatures. These results 
provide a necessary tool for a general application as well 
as for further developments of the RH approach and the 
study of Fermi-edge singularities or related phenomena 
in nonequilibrium or nonadiabatic situations. 
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APPENDIX A: TIME-REPRESENTATION OF R 

The presented method requires that the time- 
representation of R, the matrix elements of the many- 
body operator R between single-particle states, is diag- 
onal, R = R(t). For clarity and completeness, we here 
shortly state the conditions necessary for this, summa- 
rizing the arguments of Ref. (see also Refs. y| and, in 
particular, 0). For a more extended discussion we refer 
to the cited references. 

In most cases of interest, R is connected to the evo- 
lution operators U(tf) and Uo(tf) for the Hamiltoni- 
ans H and Hq, respectively. For instance, in the case 
of the Fermi-edge singularities, R represents the over- 
lap between perturbed and unperturbed states, R — 
UQ(tf)U(tf). On the other hand, for the shot noise 
spectrum of a tunneling barrier, R can be the moment 
generating function for the distribution of charge trans- 
fered out of electrode j ' = 1 during the time t = 
to tf. If Qi measures the charge of this electrode, 

Since the Hamiltonian is quadratic, the action of U (tf) 
on the many-body states can fully be described by its 
action on the single-particle states (see also Appendix[5J| , 



u(t f K\ 



(Al) 



where cr(e, e') is a N x N matrix in channel space, and 
|) the true vacuum. The matrix a can be related to the 
single-particle scattering matrix S(t, E) for a particle of 
energy E calculated for the instantaneous value of the 
potential M(t). In general, this relation is complicated, 



but becomes simple when the following condition of adi- 
abaticity is met: 



dS- 1 dS 
dt dE 



< 1. 



(A2) 



This equation encodes nothing than the condition that 
the scattering matrix has to vary slowly during the char- 
acteristic scattering time of a particle (the Wigner delay 
time). In this case, the connection to the matrix a is 
simple, 



°j,j'( e , e ') = Sj,j'(e,e'), 



(A3) 



where 

S jtj > (e, e') = \ / dt S Stj , (t, E)e^\ (A4) 

with E = (e+e')/2 and i/j the density of states in channel 
3- 

The time representation R(t) is obtained by a similar 
Fourier transform of i? JJ (e,e') = (\a e ^Ra\, ,|). 



Rj,f{e,e') = 1 j dtR jd ,(t,E)e 



i(e-e')t 



(A5) 



As noted above, the definition of R involves a or simple 
combinations of a, so that the same condition of adia- 
baticity applies. The transformations in the RH method 
do not involve the dependence of these matrices on E. At 
low temperatures and in equilibrium, the relevant physics 
is restricted to energies close to the chemical potential /i, 
so that we can fix E — /i. For nonequilibrium situations, 
the treatment remains valid as long as the dependence of 
R on E is weak. 



APPENDIX B: DERIVATION OF EQ. © 



In this appendix, we demonstrate the equality 



Tr[pR]=det(fR +(!-/)), 



(Bl) 



stated in Eq. © that connects a many-body description 
on the left-hand side to a single-body representation on 
the right-hand side. R is the many-body operator of in- 
terest, p the initial many-body density matrix at t = 0, 
p ex e~@ H ", and Tr the trace over all many-body states. 
On the other hand, R — {Rij} is the N s x N s matrix 
between the N s possible single-particle states in the sys- 
tem, R^ — (\aiRa^\) (for simplicity of notation, we here 
label single-particle states by a single index i only - in the 
main text this is split into energy e and channel indices 
j). |) is the true vacuum, not the ground state, / the 
corresponding Fermi function, and det the determinant 
over the N s x N s matrices. 

The Hamiltonian H is quadratic in the single-particle 
operators for all times t. At time t — and at zero 
temperature, the density matrix projects on the ground 
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state |0), which is a Slater determinant in the single- 
particle states = a\\). If the many-body operator 
R conserves particle numbers, it is, therefore, fully de- 
scribed by its action on the single-particle constituents 
of the ground state, represented by the matrix R. Ex- 
plicitly, for a system of N < N s particles, the ground 
state is |0> = (iV!)- 1 /2 £,(-1)^(1)) ® • • • <8> |W(iV)>, 
where the sum runs over all permutations 7r, (— 1)^ is 
the sign of the permutation, and where we assume that 
i = 1, . . . , N labels the N single-particle states of lowest 
energy. We have R\ipi) = Ri\i>i) = £j=i R ji\'*Pj)> where 
Ri is the restriction of R to the single-particle state i. 
Further R(\i/ji} <8> \ip 2 ) ® • ••) = RMi) ® #2^2) O ■ ■ • , 
from which follows that 



(0|fl|0) 



1 



E(- 



1 



(^V(l)l ® (^TT'(iV)l 



X ^[l^Tr(l)) ® IV'TT(iV)) 



<8> 



Jl,...,JJV=l 



[ 

Sis'- 1 '*" 

7T7T 7 

( — l) 7r i2 l7I -( 1 ) 



E 



'(I)tt(I) • 
Rntt(N)- 



R 



TV'(N)7T(N) 



(B2) 



This is the determinant of i?, taken over the occupied 
states only, and thus corresponds to the right-hand side 
of Eq. JEJ|. 

The generalization to nonzero temperatures is imme- 
diate: Every excited state of TV particles yields a similar 
determinant of R for a different selection of N particles 
tpj. The density matrix p weights these states with the 
Boltzmann factors oc e~P H °. Since these are exponen- 
tials of single-particle operators, they just provide the 
weights e _/3cj for the occupied single-particle states that 
sum up to the Fermi function when summing over all ex- 
cited states. Therefore, Eq. IjBlfl remains valid at nonzero 
temperatures. 
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FIG. 1: Integration contour in the f-plane. The wiggled line 
indicates the branch cut for Y(t), the isolated dots the poles 
of the integrand. The limit r — * oo is assumed. 




FIG. 2: Regions s = 1, ... ,4 of definition of the functions 
Wj (Eq. 133H ). Branch cuts are indicated by the wiggled 
lines, (a) Zero temperature case, (b) Modification for finite 
temperatures with additional branch cuts at Imt = ±/3/2. 



